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Coherent states and their generalisations, displaced Fock states, are of fundamental importance 
to quantum optics. Here we present a direct observation of a classical analogue for the emergence 
of these states from the eigenstates of the harmonic oscillator. To this end, the light propagation 
in a Glauber-Fock waveguide lattice serves as equivalent for the displacement of Fock states in 
phase space. Theoretical calculations and analogue classical experiments show that the square-root 
distribution of the coupling parameter in such lattices supports a new family of intriguing quantum 
correlations not encountered in uniform arrays. Due to the broken shift-invariance of the lattice, 
these correlations strongly depend on the transverse position. Consequently, quantum random walks 
with this extra degree of freedom may be realised in Glauber-Fock lattices. 
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Since their introduction by Glauber in 1963, coher- 
ent states have been the subject of extensive research 
within the framework of quantum optics [T] . The average 
positions and momenta of these minimum-uncertainty 
wavepackets are known to follow the motion of a clas- 
sical oscillator, thereby establishing an important bridge 
between quantum and classical mechanics Coherent 
states arise either as eigenkets of the annihilation oper- 
ator or from a displacement of the ground state of the 
quantised harmonic oscillator in phase space IJ. In gen- 
eral, if displacements of the oscillator eigenstates (termed 
Fock states or number states) are considered, a more gen- 
eral class of states, so-called displaced Fock states (DFS), 
can be obtained [31 13] . These states are of great relevance 
to many areas of quantum optics, the probably most 
important one being the direct measurement of Wigner 
functions [5], which has been successfully performed on 
propagating coherent states [5] , on single photons in cav- 
ities T and on motional states of trapped atoms [5]. 
Furthermore, DFS constitute the eigenstates of Jaynes- 
Cummings systems with coherently driven atoms [9] , and 
recently, entangled DFS have been proposed for quantum 
dense coding [TU]. DFS have been successfully generated 
by superposing a Fock state with a coherent state on a 
beam splitter [11 . However, due to the difficulties in gen- 
erating pure Fock states of higher orders, this approach 
is limited to the lowest-order DFS. To our knowledge, a 
direct observation of the genesis of these states has also 
not been possible to date. 

Quite recently, an optical system has been proposed 
which allows for a direct observation of a classical ana- 
logue for the displacement of Fock states |12] : A photonic 
lattice of evanescently coupled waveguides [13], with a 
square-root distribution of the coupling between adjacent 
guides. In these Glauber-Fock photonic lattices, every 
excited waveguide represents a Fock state and the spatial 



evolution of the light field corresponds to the probability 
amplitudes of the DFS in the number basis. Thereby, the 
emergence of these fundamental states and the underly- 
ing displacement process can be directly visualised. As 
no collapse of the wavefunction occurs for classical light, 
the displacement can be observed for a wide range of dis- 
placement amplitudes simultaneously. 
In this letter, we present the first experimental realisa- 
tion of a Glauber-Fock photonic lattice and directly ob- 
serve the classical analogue of Fock state displacements 
up to the fourth oscillator eigenstate. This arrangement 
is implemented by direct femtosecond (fs) laser waveg- 
uide inscription in fused silica [14j , whereby the required 
coupling distribution is achieved by a controlled variation 
of the distance between neighbouring waveguides. 
From another perspective, Glauber-Fock lattices provide 
a fertile ground for quantum random walks (QRWs) of 
correlated particles, as we show in the second part of this 
letter. If pairs of identical non-interacting photons are 
launched into a uniform waveguide lattice, the quantum 
interference of all possible paths leads to correlation of 
the photons [15) , thereby enabling continuous-time cor- 
related QRWs in a state space much larger than for a sin- 
gle photon [1^ . The correlation of identical particles has 
also been analysed for disordered lattices exhibiting An- 
derson localisation [T7] and for Bloch oscillations occur- 
ring in lattices with a linear gradient in their propagation 
constant [18] . As all these lattices are shift- invariant and 
infinite in the transverse dimension, the possible trajec- 
tories in a QRW are independent of their starting point. 
Breaking this invariance and introducing a boundary can 
thus embed the additional degree of freedom of transverse 
position into the QRW. We therefore analyse the correla- 
tions for pairs of separable and path-entangled photons as 
well as of fermions in a semi-infinite Glauber-Fock lattice 
and show how their correlation patterns uniquely depend 
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FIG. 1: (Color Online) Schematic view of a Glauber- Fock 
lattice of A'^ waveguides. 



on the input position. 

In order to observe the displacement of Fock states in 
the optical domain, one requires a lattice of single-mode 
waveguides, whose coupling coefficients between adjacent 
elements vary with the square-root of the site labelling 
index n 1121: 



dz 



+ Cn(j}n-1 + C„+10n+l = 0; C„ = ^^Cl . (1) 



Thereby, 4>a denotes the modal field amplitude in guide 
n, z is the longitudinal coordinate and Ci is the cou- 
pling strength between the first two guides (see Fig. [T]). 
We represent the oscillator eigenstates by single waveg- 
uides, i.e., the Fock state |k) shall correspond to the sit- 
uation when only guide k is excited: = <5nk- If light 
is launched into this single site at z = 0, the light prop- 
agation along z will map the displacement of the Fock 
state |k) along the imaginary axis of the quadrature phase 
space. More specifically, the field amplitudes evolve anal- 
ogous to the matrix elements of the unitary displacement 
operator D{a) = exp{aa^ — a*a) [Tl [T^: 



(/.n(z) = (n|I)(zCiz)|k), 



(2) 



with 

a|n) 



7,(t) 



being the ladder operators of the oscillator: 



1), 



|n) — y/n+ l\n + 1). Hence, after 
the light intensity distribution 



propagating a distance 
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D{iCiz) \k) is equivalent to the number 

distribution of the DFS D{iCiz) |k) [7. 

In the weak coupling regime, C„ depends exponentially 
on the distance between the guides: Cn = Ci cxp[— (dn — 
di)/ with K and di being fit parameters, if Ci is prede- 
termined in] . Consequently, the coupling dependence of 
Eq. ([T]) is readily achieved by inscribing the waveguides 
with = di — K\og{y/n) as the distribution of separa- 
tions. Measuring the dependence of coupling on waveg- 
uide separation at a wavelength of A = 633nm, we found 
the parameters di = 23/im and k, — 5.5^m for a desired 
coupling of Ci = 0.37cm^^. Using these results, we in- 
scribed a Glauber-Fock lattice with = 59 waveguides 
and 10cm length in fused silica, corresponding to a max- 
imum displacement amplitude a = 3.7i. We employed 
fluorescence microscopy to directly observe the intensity 
evolution of the injected light [19J and imaged the out- 
put intensity patterns onto a CCD camera. Figs. [2ja- 
d) present numerical results obtained from integrating 
Eq. ([!) for several different input sites while the experi- 
mental data is shown in Figs. [2]je-h), matching the theo- 
retical expectations very closely. The results clearly map 



the genesis of coherent states with their typical Poisson 
distribution from the displacement of the ground state of 
the oscillator (a, e) as well as the generation of higher- 
order DFS from its higher eigenstates (b-d, f-h). The 
intensity distribution features k -I- 1 maxima when site k 
is excited, in agreement to the characteristic oscillation 
with k -t- 1 maxima in the number distribution of the DFS 
I){a) \k). Interestingly, these oscillations are analogous 
to the oscillations of the Franck-Condon factors, govern- 
ing the transition probabilities between vibrational lev- 
els during an electronic transition in molecules [H [50]. 
In that context, the DFS correspond to the vibrational 
levels of the excited state, while the levels of the elec- 
tronic ground state are undisplaced oscillator eigenstates. 
Note, that the right-hand boundary at n = N — 1 is 
not reached by the propagating light. Thus, the lattice 
can be considered as effectively semi-infinite, correspond- 
ing to the semi-infinite set of Fock states. The largest 
coupling reported for laser written waveguide arrays is 
Cn-1 (X 5.5cm~"'^ |21j . which limits the maximmn num- 
ber of waveguides for the given Ci to TV c>c 200. 
As these measurements demonstrate, Glauber-Fock lat- 
tices support an optical emulation of DFS for a wide 
range of parameters, providing direct insight into their 
generation characteristics. The success of such an optical 
emulation clearly highlights the wave nature of DFS. As 
governed by Eq. ([2|, the phases of the Fock coefhcients 
are also encoded in the modal amplitudes. Hence, a full 
reconstruction of the DFS could be achieved by interfer- 
ometric phase retrieval at the end of the lattice, e.g., by 
superposing the output field with a reference wave and 
phase-stepping [52] . 

Let us now turn to QRWs of correlated photons in 
such Glauber-Fock lattices. As evident from Fig. [2] the 
intensity distribution of a classical light field is highly 
characteristic for its input waveguide. Hence, the same 
feature applies to the output probability distribution of 
a single photon. Figs. |3][a-c) show the expected distribu- 
tion for a single photon excitation of the first three sites. 
If two indistinguishable photons propagate in the lattice, 
all possible paths will interfere, and one can therefore ex- 
pect correlation patterns which are unique for each com- 
bination of input positions, in contrast to the correlations 
arising in uniform lattices where all waveguides are em- 
bedded in identical coupling environments [151 116) . 
At first, we consider input states of separable photons, 
where two photons are launched into the waveguides k 
and 1: |^'o) — a\jx\ |0). Here, a^'^ denotes the boson an- 
nihilation (creation) operator for a photon in guide k and 
|0) is the vacuum state. The probability of coincident de- 
tection of photons in guides q and r is determined by the 
photon number correlation |15| : 



Fn 



(3) 



with [/mn = )mn coupliug matrix being 

Cmn = Ci (Vn -I- l(5inn+i + V^^mn-i) for the Glaubcr- 
Fock lattice. We calculated F for a 10cm long lattice 
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FIG. 2: (Color Online) Light propagation in a Glauber-Fock lattice, (a-d) Calculated intensity evolution and output intensity 
profiles for the input sites (a) k = 0, (b) k = 1, (c) k = 2 and (d) k = 4 representing the displacement of the Fock states |0), |1), 
1 2) and |4), respectively, (e-h) Experimental fluorescent images of the intensity evolution and nearfield images of the output 
facet for a single-waveguide excitation of these sites with A = 633nm. All images have been normalised to their respective peak 
intensity. For generality, the theoretical images are shown for equidistant waveguides, while in the experiments the sites are 
distributed according to dn ~ di — Klog(y^). 
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FIG. 3: (Color Online) Photon correlations in a Glauber-Fock lattice, (a-c) Calculated output intensities for the input sites 
(a) k = 0, (b) k = 1 and (c) k = 2. (d-f) Calculated photon correlation for the input |*o) = ala} |0) with (d) (k, 1) = (0, 1), 
(e) (k, 1) = (1,2) and (f) (k, 1) — (0,2). (g-i) Measured classical output intensities at A = SOOnm. (k-m) Measured classical 
estimates. All images have been normalised to their respective peak values. 



with Ci = 0.36cm^^ and a set of input configurations 
involving the first three sites [Figs, [sjjd-f)]. The corre- 
lation patterns, and thereby the trajectories of a cor- 
related QRW, show typical bosonic bunching behaviour: 
on-diagonal peaks, corresponding to a high probability of 
detecting both photons in the same output region. Un- 
like in uniform arrays [T3] , the correlations are unique for 
each input state and depend critically on its transverse 
position, not only on the interleaving distance among the 
excitation sites. This is due to the extra phase the pho- 
tons acquire during reflection at the boundary as well as 
due to the specific Glauber-Fock coupling distribution. 



The correlation function in any photonic lattice can be es- 
timated with classical light by intensity correlation mea- 
surements |23| . We performed such measurements with 
coherent light beams in a second lattice of 59 waveg- 
uides, designed for A = SOOnm. In this spectral range 
the fabrication parameters were found to be k = 10.7/xm 
and c?i = 34/.tm for Ci = 0.36cm^^. Figs. [3]^g-i) show 
the measured output intensities from single-waveguide 
excitation, clearly exhibiting the typical number distri- 
bution of DFS. Evidently, they closely match the theoret- 
ical data presented in Figs. |3j a-c), apart from a slightly 
increased width of the light distribution, which can be 
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FIG. 4: (Color Online) Expected correlation patterns of path- 
entangled photons and separable fermions. (a-c) Photon cor- 
relation for the symmetric boson input state j ) . (d-f ) Same 
for the antisymmetric boson state I^&cT)- (g-i) Correlation 
function of two fermions being launched into sites k and 1. 



attributed to the growth of the coupling being slightly 
stronger than the desired -yn-dependence. The classi- 
cal estimates T^'^^ are obtained by launching two coher- 
ent laser beams of equal amplitude and adjustable rela- 
tive phase into the sites k and 1, measuring the intensity 
correlation among the waveguides at the output, aver- 
aging over 60 random relative phases and deducting in- 
tensity products from single-waveguide excitation. For 
a detained description of the method, we refer to [23] 
(see Eqs. (7) and (8) therein). The results obtained by 
this scheme are presented in Figs, [sf^k-m). Despite a sig- 
nificant level of noise, most likely due to fluctuations of 
phase and amplitude of the classical light fields, the es- 
sential features of F are well recovered by the classical 
estimate, namely the positions of its maxima and the 
"valleys" of destructive quantum interference. The shift 
towards larger waveguide indices is due to the increased 
spread of the light distribution as observed for single- 
waveguide excitation [Figs, [sjg-i)]. 

Finally, we theoretically study input states of path- 
entangled photon pairs, where both photons are sent 
into either of the guides k and 1 (k ^ 1): |^o) ~ 
^(a|j^ ± a^) |0). These biphoton states are the lowest- 
order A^OOA^-states, a class of maximally entangled A^- 
photon states 24 . Upon propagation of such a state 
in a photonic lattice the photon correlation yields pS] : 
Fqj, = |C/qk^^rk ± C^qiC^ril • The corrclatious have been 



calculated for the lowest order Fock states for the sym- 
metric (-I-) as well as for the antisymmetric state (— ) 
[Figs. |4]Ja-f)]. As before, the correlation maps are highly 
distinct for each particular input configuration and be- 
come more involved when higher excitation sites are used. 
To gain a more thorough insight into the nature of these 
correlations, it is worthwile to compare them to the corre- 
lation of product states of identical bosons [Eq. ^] and 
fermions in Glauber-Fock lattices. As recently demon- 
strated, fermionic correlations can be accurately simu- 
lated in birefringent waveguide lattices, provided that the 
two polarization states are entangled across different sites 
Introducing the annihilation(creation) operators 

and considering initial 

') (k 7^ 1) , the fermionic 



bn^ for fermions in a lattice Ell 



states of the type 
correlation 



function 

2 



/ 

yields [17: Pg^ = (fe^6j^6q) = 
UqkUii — C/qiC/rkT- The results displayed in Figs. k[g-i) 
show a strong anticorrelation: Off-diagonal peaks domi- 
nate the correlation map and there is zero probability to 
find both particles in the same channel, as required by 
Pauli's exclusion principle. Quite remarkably, the cor- 
relation patterns for the entangled photons [Figs. Qa- 
f)] reveal a composition of bosonic bunching [compare 
Figs, [sj^ d-f)] and fermionic antibunching features, similar 
to the situation in disordered lattices [T7] . In some cases 
the main peaks follow a bosonic behaviour(c, d, e), while 
in others the fermionic anticorrelation prevails (a, b, f). 
However, as the entangled photons are still bosons, they 
never exhibit the strong on-diagonal zero trench char- 
acteristic for fermions. Note that a classical estimate 
for the correlation of path-entangled photon pairs can be 
obtained from intensity correlations of two coherent light 
beams as well, but requires a precisely controlled relative 
phase [25]. 

In conclusion, we directly observed a classical analogue 
for the displacement of Fock states by monitoring the 
propagation of classical light waves in Glauber-Fock pho- 
tonic lattices. This demonstrates that essentially wave 
mechanics governs the displacement process. As every 
waveguide is subjected to a different coupling environ- 
ment, correlations of identical particles evolving in such 
a lattice are highly characteristic for each input configu- 
ration. Therefore, quantum random walks exploiting the 
additional layer of complexity associated with this degree 
of freedom seem in reach. 
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